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The nucleation and structure of silicon nanocrystals formed by different preparation conditions 
and silicon concentration (28 - 70 area %) have been studied using Transmission Electron Mi- 
croscopy (TEM), Energy Filtered TEM (EFTEM) and Secondary Ion Mass Spectroscopy (SIMS). 
The nanocrystals were formed after heat treatment at high temperature of a sputtered 10 nm 
thick silicon rich oxide on 3 nm Si02 layer made by Rapid Thermal Oxidation (RTO) of silicon. 
Nanocrystals precipitate when the excess silicon concentration exceeds 50 area %. Below this per- 
centage amorphous silicon nanoclusters were found. In-situ heat treatment of the samples in the 
TEM showed that the crystallization requires a temperature above 800°C. The nanocrystals pre- 
cipitate in a 4 nm band, 5 nm from the Si substrate and 4 nm from the Si02 sample surface. The 
silicon nucleates where the excess Si concentration is the highest. The top surface has less excess 
Si due to reaction with oxygen from the ambient during annealing. The Si02-RTO layer is more Si 
rich due to Si diffusion from the Si02-Si layer into RTO. Twinning and stacking faults were found in 
nanocrystals with 4-10 nm in diameter. These types of defects may have large effects upon the us- 
ability of the material in electronic devices. Both single and double twin boundaries have been found 
in the nanocrystals by high resolution transmission electron microscopy (HRTEM). Image simula- 
tions were carried out in order to obtain more information about the defects and nanocrystals. The 
stacking faults are extrinsic and located in the twin boundaries. 

PACS numbers: 



I. INTRODUCTION 

Nanoscaled electronic devices have attracted much 
attention, especially related to MOS (Metal-Oxide- 
Semiconductor) devices used for memory storage 
applications 12 . The improved electrical properties and 
the possibility to scale electronic devices down consid- 
erably are important and useful. Silicon nanocrystals 
(NCs) embedded in silica can potentially be used for 
various applications 3 ' such as nanocrystal mem ory cells 4 , 
photon converters, optical amplifiers etcP^. It is as- 
sumed that the replacement of the bulk floating gate with 
nanocrystals will in turn result in longer retention, lower 
gate voltage and lower power consumption^. The sepa- 
ration between the nanocrystals can prevent charge loss 
laterally and results in short writing times at lower volt- 
ages and improved reliability^. Accurate control of the 
array of nanocrystals is very important, since even 1 nm 
change in tunnel distance can affect the write and erase 
time^. Both injection and retention of electrons in these 
devices are very sensitive to the size, distribution, inter- 
faces and electronic structure of the nanocrystals. It is 
desirable to make nanocrystals that are less than 10 nm 
in diameter so that the Coulomb blockade effect becomes 
prominent at room temperature^. When the dimen- 
sions approach the atomic scale, significant changes occur 



in the electronic, optical and thermodynamic properties 
compared to bulk materials^. Therefore it is important 
to get information about the nucleation of amorphous 
and crystalline nanoclusters, to characterize their atomic 
and electronic structure, the nanocrystal interfaces, and 
the defects within and around the nanocrystals. 

Twins are among the most common defects in SP 
and depend on the stacking fault energy of the mate- 
rial, the surface stresses and surface orientations^. Low 
stacking fault energy results in an increase in twinning^. 
Nanocrystal twinning can be due to stresses in and on the 
surface of nanocrystals, coalescence of smaller nanocrys- 
tals or it can occur during growth and heat treatment. 
Wang et alP^ studied nanocrystals made by ion implan- 
tation of Si in a 1/im Si02 film. They found twinning and 
stacking faults in 90 % of the nanocrystals with a crystal 
size larger than 5 nm. They claimed that nanocrystals 
smaller t han 5 nm will not contain any twinning or stack- 
ing faults^M!. Perrey et alP studied silicon nanocrystal 
defects, precipitated from hydrogenated amorphous thin 
Si films (a/nc-Si:H). Both twinning and stacking faults 
were found in the silicon nanocrystals with sizes 1.5 - 
5 nm in diameter. Few similar studies have been done 
on nanocrystals made by sputtering in ultra-thin oxides 
(less than 20 nm in thickness). 

Most previous studies of 3-12 nm nanocrystals have 
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been carried out on thick silicon oxides (500nm - 1/zm) 
with implantation doses bet ween 2*10 16 - 3*10 17 Si + 
cm- 2 (3-50 at. %) ii | im4 | istf I5] Si nanocrystal forma- 
tion, size and distribution in thin oxides have previously 
not been studied in detail. The nucleation of Si nanocrys- 
tals depends on the density of the SiC>2 layer and defects. 
This varies with different sample preparation techniques. 
Nanocrystals have previously been made by methods like 
sputtering, sol-gel, ion implantation and other chemical 
processes^. The nanocrystals in the present work were 
intended to be formed in a 10-15 run thin oxide and close 
to the surface. The concentration and nanocrystal size 
in ion implanted samples ha ve pre viously been shown to 
have a Gaussian distributio n 1 14 * 18 ^ with a peak position 
further down in the oxide. Ion implantation could not be 
used on such thin oxides, due to less control of the aver- 
age cluster size and density^. Therefore sputtering was 
used. The nucleation of nanocrystals, distribution, crys- 
tal size and the defects were studied in the present work 
by High Resolution Transmission Electron Microscopy 
(HRTEM), Energy Filtered TEM (EFTEM) and Sec- 
ondary Ion Mass Spectroscopy (SIMS). In order to get 
detailed information about the nanocrystals and their de- 
fects, experimental through focus series were compared 
to simulated images. 



II. EXPERIMENTAL 

The samples were made by growing a 3 nm thin layer 
of SiC>2 on a p-type silicon substrate by Rapid Thermal 
Oxidation (RTO) at 1000°C for 6 sec. Then a 10 nm 
layer of silicon rich oxide was sputtered from a Si0 2 :Si 
composite target onto the RTO and heat-treated in a N 2 
atmosphere at 1000 - 1100°C for 30-60 min. Different 
area percentages (area coverage) of Si:Si02 (6, 8, 17, 28, 
42, 50, 60, 70 area % corresponding to 4, 5, 11, 18, 33, 
40, 46 at. %) were used to produce different silicon su- 
persaturations in the oxide. The samples dicussed in this 
paper are given in Table [T] Cross-sectional TEM samples 
were prepared by ion-milling using a Gatan precision ion- 
polishing system operated at 5 kV gun voltage. Silicon 
nanocrystals were observed by HRTEM and Energy Fil- 
tered imaging of the plasmon peak. The HRTEM and 
Energy Filtered TEM (EFTEM) were performed with 
a 300 keV JEOL 3100FEF microscope with an Omega 
energy filter. The spherical (C s ) and chromatic aberra- 
tion (C c ) coefficients of the objective lens were 0.6 mm 
and 1.1 mm, respectively. The point to point resolution 
was 0.174 nm at Scherzer focus (-37 nm), and the mini- 
mum probe diameter was 0.2 nm. The energy resolution 
at 300 keV was 0.78 eV, and was previously determined 
experimentally^!. The energy dispersion of the Omega- 
filter was 0.85 /im/eV at 300 keV. Through focus series of 
HRTEM images of crystal orientation [110] were recorded 
with an objective aperture large enough to include the 
220 reflection (0.192nm) corresponding to 5.6 nm -1 (11 
mrad) radius. Ten images with a 10 nm difference in 



focus were obtained around Scherzer condition. The im- 
ages were compared to simulated focus series made by 
using the MacTempas computer program. In-situ heat- 
ing of sample 70asd (see Table [T]) was in addition carried 
out with the JEOL 3100FEF microscope using a heat- 
ing holder. EFTEM- Spectral Imaging (EFTEM-SI) was 
acquired by scanning from 2 eV to 30 eV energy losses, 
with an energy slit of 2 eV. Electron energy loss spectra 
were aquired with a 197 keV JEOL 2010F microscope 
equipped with a Gatan imaging filter and detector. Sec- 
ondary ion mass spectrometry (SIMS), using a Cameca 
IMS 7f with a 25 nA Oj, 0.5 keV primary beam was 
used to obtain depth profiles of Si secondary intensities. 
The ion beam was scanned over an area of 200 x 200 
/im 2 . Positive secondary ions were collected from a cir- 
cular area with a diameter of 62 fim from the center of 
the crater and counted with an electron multiplier. Af- 
ter the SIMS measurements the depth of the craters were 
measured with a Dektak 8 stylus profilomctcr, for con- 
verting the sputter time to depth. The depth calibration 
was made on the assumption of a constant erosion rate. 

TABLE I: The samples discussed in this paper, presented with 
Si concentration, heating time and temperature. 
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III. RESULTS AND DISCUSSION 

A. Nucleation of crystalline and amorphous 
nanoclusters 

The nucleation of nanocrystals with different preparation 
conditions was studied by HRTEM and EFTEM. No sil- 
icon nanocrystals were found in sample 28 (see Table [i]). 
At silicon fraction of 42 % (sample 42) 4-6 nm bright, 
isolated areas with abrupt interfaces are visible in the 
EFTEM image shown in Figure [I] The graph in the left 
corner shows how the contrast changes across one of the 
bright area in the middle of the image. The bright areas 
are located about 5 nm from the silicon substrate and 4 
nm from the specimen surface (Si02/glue). The EFTEM 
image was produced by using the plasmon Si peak at 16.8 
eV. HRTEM images show no visible nanoclusters or any 
contrast differences in this area. The bright, isolated ar- 
eas in Figure [I] of sample 42 are probably amorphous Si 
nanoclusters, because the local changes in image inten- 
sity cannot be explained by variations in sample thick- 



3 



ulk-Si 



Si02-Si\ | 




FIG. 1: EFTEM image (at 16.8 eV) of amorphous nanoclus- 
ters in sample 42. The arrows indicates the amorphous nan- 
oclusers which are Si rich. The graph in the left corner shows 
the differences in intensity over the bright dot in the middle. 



ness. Since no nanoclusters were found in sample 28, the 
excess silicon had either oxidized, precipitated to small 
amorphous nanoclusters or are still silicon atoms or ions 
in solid solution. EFTEM images of the surface plasmon 
peak of Si were also acquired, but no such bright areas 
were visible. 

HRTEM and EFTEM images of samples 50 and 70 
show areas with lattice fringes indicating crystalline nan- 
oclusters (see Figure [2). The nanocrystals in these 
samples precipitate at the same distance from the sub- 
strate as the amorphous nanoclusters in sample 42. The 
nanocrystals first grow spherically up to about 4 nm in 
diameter, following an elongated growth laterally in the 
<111> directions (see Figure [2]). The nanocrystals have 
the same diamond type structure as bulk silicon with lat- 
tice parameter 0.54 nm. These observations show that 
increasing silicon concentration, annealing time or tem- 
perature induce nanocrystal growth. 

Sample 70asd was examined with HRTEM and 
EFTEM-SI. No lattice fringes were observed in the 
HRTEM images. With EFTEM-SI small bright areas 
were observed. After an in-situ heat treatment in the 
TEM for 1 hour at 800° C, still no indications of crys- 
talline nanoclusters were found (see Figure[3| . The bright 
areas in the EFTEM images are probably small amor- 
phous Si nanoclusters. During heat treatment, the Si 
atoms in the nanoclusters diffuse and form larger areas of 
amorphous Si. Crystalline nanoclusters were only formed 
after heating at a temperature above 800° C, see Figure 

m 

EDS spectra of the Si0 2 -RTO and the Si0 2 -Si layer in 
sample 70 are presented in Figure |4^. The location of the 
EDS measurements are shown in Figure [2^,. Both EDS 



8) Si02 Si02 
top layer NCs (RTO) 

(Si02-Si) 



1 




FIG. 2: Images of sample 70. a) HRTEM image, b) EFTEM 
image, with the plasmon peak of silicon (16.8 eV). The num- 
bers indicate the area where the EDS and EELS spectra were 
obtained (see Figure j4|). 



spectra were acquired with the same spot size and acqui- 
sition time. The Si-Ka EDS peak have almost the same 
intensity for both regions, but slightly higher from the 
Si0 2 -Si area than the RTO. The latter area (RTO) has 
slightly higher sample thickness compared to the Si02-Si 
(area 2) layer. Therefore the differences in EDS inten- 
sities indicate that the RTO layer is probably not only 
composed by Si0 2 , but contain high amounts of pure Si. 

The increase in silicon concentration in the RTO com- 
pared to pure Si0 2 can also be seen from the EELS spec- 
tra in Figure |4Jd. These spectra were acquired from the 
Si0 2 -RTO (area 1) and the top layer (area 3), see Figure 
[2]}. The RTO is slightly thicker than the top layer, and 
therefore has a higher plasmon peak (as seen in Figure 
[4}d). There is a small difference (1.2 eV) in plasmon peak 
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FIG. 3: EFTEM-SI of sample 70asd after an in-situ heat 
treatment for 1 hour at 800° C. The image is made with the 
plasmon peak of Si (16.8 eV). 
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energy between the RTO (22.8 eV) and the top layer 
(24.0 eV). The difference is within the standard devia- 
tion of the spectrometer (1.5 eV), but the experimental 
observed energy shift may indicate trends. The plasmon 
peak energy of pure silicon is 16.8 eV, while pure Si02 
has 24 eV. A shift towards lower plasmon peak energy 
may therefore indicate an increase in silicon concentra- 
tion in the oxide. The plasmon peak of the RTO exhibits 
in addition a broader width in the range of 19-25 eV. 
This is due to a significant increase in the pure silicon 
contribution (at 16.8 eV) in the oxide^. 

SIMS data in Figure [4]: shows the Si secondary inten- 
sity profile of sample 70 and 70asd. The Si concentration 
maximum is located 7 nm from the sample surface. This 
could be due to oxygen diffusion from the surface into 
the oxide. There is also an increased Si concentration in 
the Si02-RTO (area 1) than may be expected from pure 
Si02- This could be due to silicon diffusion from the 
Si0 2 -Si layer (area 2) into the Si0 2 -RTO layer. Sam- 
ple 70 and 70asd have different Si secondary ion profiles. 
This suggests that diffusion has occurred both before and 
after heat treatment. The observed Si secondary inten- 
sity is higher in the oxide than in the bulk due to a ma- 
trix effecP^l. The amorphous and crystalline nanoclus- 
ters both nucleate at the concentration maximum. This 
maximum is located about 7 nm below the surface and 
is due to both oxygen and silicon diffusion before and 
after the heat treatment. When sample 70 and 70asd is 
exposed to air, oxygen diffuses into the top layer (area 3) 
and additional silicon diffuses from the Si02-Si or from 
the Si-substrate into the RTO. This is in agreement with 
both the EDS, EELS and SIMS data in Figure [4] 
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FIG. 4: EDS, EELS and SIMS measurements on sample 70 
and 70asd. The location of the measurements are shown in 
Figure [2^,. a) EDS spectra of the silicon Ka line are from 
the Si0 2 -RTO (area 1) and the Si0 2 -Si (area 2). b) EELS 
spectra of the Si02-RTO (area 1) and the top layer (area 
3). The two spectra are normalizes with regard to the zero 
loss peak, c) SIMS data from sample 70 and 70 asd. The 
spectra are normalized so that the Si matrix has a secondary 
ion intensity of 1.0. 



B. Analysis of HRTEM images 

In HRTEM images of thin silicon specimens, bright dots 
at Scherzer conditions correspond to the tunnels in the 
structure^. In thicker regions the bright dots in the im- 
ages correspond to the atomic columns. To distinguish 
between these two cases in a perfect crystal, is impossible 
based only on the image contraslP^. Figure [5^, shows a 



FIG. 5: a) A (111) twin boundary in a diamond type struc- 
ture, projected along [Oil] and b) a simplified projection seen 
at limited resolution. 



schematic drawing of a (111) twin boundary in a diamond 
type structure in the [Oil] projection and Figure]^ shows 
a schematic simplified image at limited resolution. Since 
the twin boundary exhibits a mirror symmetry across a 
plane of tunnels (rather than atomic planes), the posi- 
tion of the mirror line in experimental HRTEM images 
may be used to determine whether the bright dots in 
an experimental image correspond to tunnels or atomic 
columnsP. The experimental HRTEM image in Figure 
[6] shows mirror symmetry across a line located between 
two rows of bright dots. Therefore the dots correspond 
to the atomic columns. The mirror symmetry across the 
two twin boundaries in the nanocrystal image in Figure 
[7] is present along a line with bright dots as indicated 
by arrows. Therefore the bright dots correspond to the 
tunnels in the silicon structure. In order to extract more 
information about the atomic structure of the nanocrys- 
tals a series of ten experimental HRTEM images were 
taken with a Af=10 nm difference in focus. The focus 
of the recorded HRTEM images was determined exper- 
imentally. This was done by optical diffractograms by 
Fast Fourier Transformation (FFT), see Figure H The 
method usi ng op tical diffractograms was developed by 
Thon in 1971* 23 * 25 -l The experimental diffractograms from 
sample 70 showed three bright rings. The radial inten- 
sity of the ring pattern is approximately proportional to 
sin 2 x(u) where x(u) is the transfer function given by^ 

x( U ) = (^AA/ u 2 ) + (|)(A 3 C sM 4 ) (1) 

A is the electron wave length, u is the distance to the 
maxima and minima in the optical diffractogram. C s is 
the spherical aberration coefficient, Af is the defocus of 
the objective lens. The maxima or minima in the diffrac- 
tograms correspond to sin 2 x(u)= 1 or 0. Equation [l] can 
be transformed to 

(n/u 2 ) = (2A)A/+(C S AV (2) 

where n odd for maxima and n even for minima. A plot 
of y=n/u 2 as function of x=u 2 gives a line that intersect 
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FIG. 6: Single-twin in a nanocrystal in sample 70. 




FIG. 7: Double-twin in a nanocrystal in sample 70. 



the y-axis at b=2AAf and with a slope a=C s A 3 . For the 
experimental HRTEM images recorded in this work the 
electron wave length A was 0.00224nm as calculated for 
U = 300 kV and the instrumental value of the spherical 
aberration coefficient C s = 0.6 mm was used. This gave a 
slope of a= 0.0067nm 4 . The intersection with the y-axis 
b=2AAf was used to determine the defocus. 

The function n/u 2 versus u 2 as given by Equation [2] 
was plotted with n=-l to n=-10. The curves are pre- 
sented in Figure [8] The square of the reciprocal distances 
(u) from the central beam to a maxima and minima of 
sin 2 x(u) in the optical diffractogram are plotted along 
the x-axis. The intersections between these points for 
a particular HRTEM image and the theoretical curves 
(with different n values), should satisfy the linear equa- 
tion [2j The two straight lines shown in Figure [8] were 
calculated with regard to HRTEM image 6 and 10 in 
Figure [9) The defocus (Af) was then determined by the 
intersection of the straight line with the y-axis b=2AAf. 
This procedure was carried out for each HRTEM image 
in the focus series, and was found to cover a variation in 
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FIG. 8: A plot of ii/m 2 as a function of u 2 . The plot is used 
to determine the focus of the HRTEM images in Figure [9] 
The diffractogram in the bottom right corner corresponds to 
HRTEM image no. 10 (see Figure^. The spots in the diffrac- 
togram from the Si substrate was used as internal standard. 
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FIG. 9: Experimental through focus series of the nanocrystal 
with a double twin structure as in Figure [7J together with two 
series of simulated HRTEM images at different thickness and 
focus. 



Af from -20nm to -llOnm. Unfortunately the measure- 
ments of the maxima and minima in the diffractograms 
have a relatively large standard deviation and this results 
in an accuracy of typically ±10nm in Af. 

The experimental focus series were compared to sim- 
ulated HRTEM images for different sample thicknesses 
with focus variations from Af= -20nm to Af= -llOnm. 
The simulated through focus series were performed us- 
ing the MacTempas computer program. The best agree- 
ment between the experimental and simulated images 
was found for thicknesses to be between 20 and 30 nm, 
see Figure [9j This figure presents two series of simu- 
lated HRTEM images as a function of focus, calculated 



for sample thickness of 20 and 30 nm. These images 
were then used to identify the relationship between the 
image contrast and atomic structure near a twin bound- 
ary. From these images it is clear that the bright dots 
in the experimental HRTEM images with a focus Af=- 
20nm to Af=-60nm correspond to the tunnels. For the 
images with other focuses (-70 to -110) the bright dots 
correspond to the atomic columns as can be deduced from 
the position of the mirror line in the image. 

The Fourier image period was calculated using equa- 
tion § for various reflections {uhki) inside the objective 
aperture^. 



A// 



(3) 



The Fourier image period is the period at which the 
first and last HRTEM images in a series are identical. 
The point resolution of the microscope was calculated 
to be X mi „=0.66Cy 4 A 3 / 4 =0.174 nm. Lattice distances 
in Si larger than the point resolution of the microscope 
are dm— 0.314nm, g?200= 0.272nm and ^220= 0.192nm. 
Therefore the objective aperture used for the experimen- 
tal images was large enough to include these 3 types 
of reflections. The Fourier image periods corresponding 
to these reflections are A//(tim)= 97.7nm, A//(u 2 oo)= 
73.8nm and A//(it22o)= 36.9nm respectively. By com- 
paring these three periods with the HRTEM images in 
Figure [9j the periods 36.9nm and 73.8nm do not satisfy 
the conditions of the Fourier image period. The first and 
last calculated HRTEM images in the two periods are 
not identical. This shows that the Fourier image period 
of the experimental image in Figure [9] is 97.7nm. 



Nanocrystal defects: Stacking faults and 
twinning 



Nanocrystals smaller than 3 nm in diameter show large 
lattice strains on planes, see Figure 10 
planes might be due to surface effects 
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nancrystal are embedded in the Si0 2 matrix. During 
nanocrystal growth, this strain might lead to both twin- 
ning and stacking faults. 

In the crystalline nanoclusters twinning and stacking 
faults occur, even in clusters down to 4 nm in size. In 
large elongated nanocrystals no defects were observed. 
Figure [6] presents a 7.5 nm nanocrystal with a single twin 
structure, as seen in the [011] projection. The twin plane 
is indicated by an arrow labeled TB. The angle between 
the atomic planes on opposite sides of the twin boundary 
is 70.5°, which agrees well with the value 70.53° between 
the {lll}-type planes in cubic materials. The projected 
shape of the nanocrystal is hexagonal. This is in agree- 
ment with the report by Wang et alP^, who also found a 
hexagonal shape of nanocrystals containing a single twin 
boundary. The nanocrystal boundary consists of {111} 



7 




FIG. 10: HRTEM image of a small (3nm) nanocrystal in sam- 
ple 50. 



and {110} crystal facets. The parallel white lines illus- 
trate a deviation in the atomic stacking characteristic for 
stacking faults. 

Figure [7] presents a 4 nm nanocrystal with a dou- 
ble twin structure and two stacking faults. The twin 
planes are {111} planes and the twin boundaries are in- 
dicated by arrows labeled TBI and TB2. Facets are 
visible along the {111} planes. The two twin bound- 
aries are parallel and only separated by 4 atomic tun- 
nels. The nanocrystal stacking faults are located in the 
twin plane and is passing through the whole nanocrys- 
tal. Both twins and stacking faults frequently occur in 
the equivalent plane, the {111} plane, in diamond type 
structures^. The atomic sequence of the stacking fault 
and twin (ABCBACBA) is shown on Figure[7] The stack- 
ing fault is extrinsic, and can be formed by inserting an 
extra atomic plane B in the perfect crystal (ABCABC 
to ABCBA). Extrinsic stacking faults can be described 
by use of shears in successive planes in a perfect crystal. 
Wang et alP^ proposed that when the stacking faults and 
twinning coexist, the projected shape becomes more ir- 
regular. The projected shape of the nanocrystals in the 
present work was not irregular. This may lead to the con- 
clusion that when stacking fault is located in the twinning 
planes, the projected shape is the same as in a stacking 
fault free crystal. 

Contrary to the samples studied by Wang et alP^ 
made by ion implantation, twinning and stacking faults 
were found in nanocrystals smaller than 5 nm. Perrey 
et alP^ studied silicon nanocrystals embedded in amor- 
phous hydrogenated silicon. They also found twinning in 
nanocrystals smaller than 5 nm in diameter. The sample 
preparation method and matrix type around the silicon 
nanocrystals seem to affect the defect concentration in 
the nanocrystals. Increased SiO^ in the matrix may im- 
ply higher surface stresses to the nanocrystals and may 
induce twinning and stacking faults. 



IV. CONCLUSION 

Nanocrystals of silicon in a thin Si02 matrix upon a 
Si substrate were studied with HRTEM, EFTEM and 
SIMS. The nucleation and their structure in silicon were 
examined in detail. The studies show that nanocrys- 
tals precipitate at a silicon fraction of 50 area % and 
above. Below this percentage amorphous silicon nan- 
oclusters were found. Both crystalline and amorphous 
nanoclusters precipitate in a 4 nm band, 5 nm from the 
silicon substrate and 4 nm from the specimen surface. 
EELS, EDS and SIMS show a higher silicon concentra- 
tion in the RTO layer and more oxygen in the top layer. 
These results may indicate that oxygen diffuses into the 
oxide close to the surface, whereas silicon diffuses from 
the Si0 2 -Si layer and into the RTO during heat treat- 
ment. The defects in the nanocrystals were studied in 
detail, and found to be twinned clusters with extrinsic 
stacking faults. The smallest nanocrystal with defects is 
4 nm in diameter. The small amorphous nanoclusters 
have a spherical shape while the crystalline nanocrys- 
tals have a more hexagonal shape with facets. After fur- 
ther heat treatment, higher temperature or silicon con- 
centration, the nanocrystals have an elongated growth in 
the <111> direction. Detailed HRTEM image simula- 
tions were carried out to get more information about the 
nanocrystals and their defects. All stacking faults found 
in the nanocrystals are extrinsic, going through the whole 
nanocrystal and is located in the twin boundary. No dis- 
locations were observed. 
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